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A  Method  for  Evaluating  a  Generalized  Fresnel  Integral 
Related  to  the  Spectra  of  Amplitude  and  Frequency  Modulated  Pulses 

A.  Statement  of  the  Problem 

Recently  questions  have  arisen  concerning  the  magnitude  of  the  inter¬ 
ference  occurring  between  adjacent  radars  at  the  present  AMRAD  site.  These 
questions  have  created  the  need  for  evaluating  the  integral 

2  . 

G(f,  O')  ^  g(t)  g-i2Trft 

-  00 

for  various  functional  forms  of  g(t)  and  for  various  ranges  of  the  magnitudes 

2 

of  the  essential  parameters  fT  and  o^T  ,  where 

t  =  time;  f  =  frequency  in  cycles  per  second;’^' 

/  2 

O'  =  the  frequency  modulation  parameter  in  radians/sec  ; 

T  =  the  pulse  duration  in  seconds;  (2) 

g(t)  =  h(t)  in  0  t  ^  T ;  h(t)  =  either  a  real  or  a  complex  function; 
g(t)  =0  in  t  <  0,  t  >  T. 

When  the  function  g(t)  is  defined  in  the  above  manner,  Eq.  (1)  can  be  re¬ 
written  as  follows; 

G(f,  a)  =  ^  h(t)  e'  dt  .  (3) 

o 


^The  variable  f  can  be  replaced,  if  the  analysis  required  it,  by  the  difference 
frequency  Af  =  f  -  f^,  where  f^  is  the  CARRIER  FREQUENCY  and  f  is  the 
frequency  at  which  the  spectrum  G(f)  is  viewed.  In  this  case,  f  is  required 
to  lie  in  a  range  corresponding  to  the  inequality  stated  in  Eq.  (6),  but  is 
otherwise  arbitrary. 
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Evaluation  of  this  integral  yields  in  general  the  spectrum  of  a  general  pulse 
h(t)  modulated  in  amplitude,  phase  and  frequency.  The  first  two  types  of 
modulation  correspond  to  variations  in  the  amplitude  and  phase,  respectively, 
of  the  complex  function  h(t).  Frequency  modulation  corresponds  to  non¬ 
vanishing  values  of  the  parameter  a.  In  the  present  paper,  only  amplitude 
and  frequency  modulation  will  be  considered,  i.  e.  ,  h(t)  will  be  assumed  to 
be  a  real  function  of  t. 

The  integral  in  Eq.  (3)  is  a  generalized  Fresnel  integral  in  the  sense 
that  for  arbitrary  functional  forms  of  h(t)  it  is  a  generalization  of  the 
integral 


(4) 


which  can  be  reduced,  with  the  aid  of  the  proper  change  of  variables,  to  a 
difference  of  two  classical  Fresnel  integrals  of  the  type  encountered  in 
Optics,  namely,  of  the  following  type: 


u 

J(u)  =  ^ 
o 


de  . 


(5) 


Now  the  behavior  of  the  special  integrals,  Eqs.  (4)  and  (5),  is  well 
known,  and  extensive  tables  are  available  for  use  in  evaluating  them.  The 
same  cannot  be  said,  however,  for  the  general  case,  when  h(t)  in  Eq  (3) 
is  other  than  a  constant.  Very  little  information  is  available  in  the  literature 
regarding  the  behaviour  of  the  generalized  integral.  In  the  absence  of  such 
information,  it  was  necessary  to  devise  an  original  approach  to  the  evaluation 
of  Eq.  (3).  The  approach  developed  by  the  author  and  employed  in  the  present 
paper  is  outlined  in  Section  C. 
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B.  Summary  of  Results 

The  integral  in  Eq.  (3)  has  been  evaluated  for  nine  different  cases  -- 

eight  of  them  specific  cases,  the  other  a  general  case.  In  each  case,  the 

2 

essential  parameters  fT  and  O'T  were  assumed  to  obey  the  following 
inequality: 


ZirfT  »  aT^  »  1  ,  (6) 

Or  equivalently,  in  terms  of  variables  A  and  A  introduced  later  in  the 
analysis,  the  following  inequality; 

A  A  1  A  _  fT  A  ^ 

^  =  TT  ^  =  ~7fr  • 

The  general  case  referred  to  above  corresponds  to  the  following 
representation  of  the  function  h(t)  in  the  interval  0  ^  t  ^  T:  h(t)  a  piece- 
wise~linear  function.  * 

In  this  case,  assuming  the  interval  0  ^  t  ^  T  to  be  divided  into  N 
equal  intervals,  each  of  length  T/N,  we  have 

h(t)  =  hjt)  in  (n-l)(l)  <  t  <  (n)(^)  ;  (8) 

where 

h  (t)  =  a  +  a,  t  .  (9) 

n' '  on  In 

The  eight  specific  pulse  shapes  which  were  treated,  together  with 
the  corresponding  definitions  of  the  function  h(t),  are  listed  below  and  also 
in  Fig.  1.  In  the  equations  below,  the  function  h(t)  is  followed  in  each  case 
by  the  associated  spectral  function  G(f,  ca)  obtained  by  evaluating  Eq.  (3) 

*A  method  for  calculating  the  FM  spectrum  of  a  pulse  approximated  by  a 
piecewipe  linear  function  is  outlined  in  Section  D. 
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with  the  aid  of  the  methods  outlined  in  Section  C.  * 

1.  Rectangular  Pulse 

h{t)  =  1  (10) 


G(f,o)  =.  i(^) 


exp  Q  [aT^  -  ZufT}] 

(f  - 


1 

T 


+  o{ 


4) 


(H) 


2.  Triangular  Pulse 

h{t)  =  '^  inOiSt:^-^ 

h{t)  =  +  T)  in  ^  <  t  <  T 


(12) 


G{f,a)  -  -{— 4- 
2tt  T 


7 


2 


exp 

(f 


[ifQ'T^/4  -  TTfT}] 

aT.2 

■  1^’ 


r  r  2  ^  'll 

exp 

i  (q-T  -  2TTfT] 

(f 

aT.2 

IT  ' 

+  0(4) 

f 


(13) 


‘•‘The  expressions  for  G(f,  O')  given  in  Eqs.  (11),  (13),  (15),  (17),  (19),  (21), 

(23),  and  (25)  were  obtained  in  each  case  by  approximating  G(f,  O')  by  the 

first  non-vanishing  term  by  the  asymptotic  expansion  of  G(f,  O')  (of  the  type 

derived  in  Section  C).  When  G(f,  Q^)  consists  of  the  sum  of  several  separate 

expansions,  taking  the  first  term  means  adding  the  first  term  (non-vanishing) 

of  each  of  the  component  expansions.  For  example,  in  Eq.  (43)  one  would 

take  the  sum  of  the  first  non- vanishing  terms  of  the  three  expansions 
MM  M 

Z  ,  Z  ,  and  Z  . 
m=o  m  =  l  m  =  2 
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3.  Parabolic  Pulse 


h(t)  =  (^)  (1  -  y) 


(14) 


G(f.a)  s  -(^) 
TT  T 


7 


exp  i  {cuT^  -  ^Trflj] 


(f 


[ifaT^  -  Z-rrfT)] 


1  exp  I  i  faT 

?  — (f 


+  0(4-)  (15) 

f 


4.  Trapezoidal  Pulse 

h(t)  =  h^(t)  =  t/r  in  0  <  t  <  t  ;  t  -  —  ^  5  ^  —  T  ; 

h(t)  =  ^^(t)  =  1  in  T_  <  t  <  T^;  =  "4~'' 


(16) 


h(t)  =  h^(t)  =  (-t+  T)/t  in  <  t  <  T 


G(f,a)  -  -( 


exp 


7 


l^i  |aT  ^  -  Z-rrfT  |J 


ar 


(f 


exp  '  2TrfT|j 


ar 


(f  - 


+  ^2 


(f  -W 


+  0(4-)  (17) 

f 


5.  Parabolinear  Pulse 


h(t)=h^(t)=ajjt  in  0<t<T;T 


T  -  e 

~Z~ 


'll 


T  .2 


h(t)  =  h,(t)  =  1  -  (— — )(t  -  -7  )  in  T  <  t  <  T  ; 


h(t)  =  62(1)  =  a^  j  (-t  +  T)  in  <  t  <  T  ; 


T  = 


T  -  € 


T  +  e 


;  a 


11  2T  -  e  -  4t  +  e 


(18) 


J 
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G(f,a)  = 


+ 


+ 


where 


6.  Sinusoidal  Pulse 
h{t)  =  sin 

G(f,c^)  ^ 

+ 

S  { 


^11 

4-n^ 


1  1 

l(— 

4ir  € 


1  exp  [ifo-T^  -  Ztt£tJ\ 

7  (£  - 

'  TT 


exp 


j^i  ^  -  Zirfr  jj 


ar 


(f  - 


TT 


exp 


-  2TrfTj] 


ar 

(f  -  -^  )^ 


+  O(^) 

f 


(19) 


/  1  X  .  ;/  11  ^  -  /  1 

(  T  Q  )  (  ■ 

TT  (2T  -  e)  e 


(tt  e)(2T  -  e) 


) 


(19a) 


(?)  =  (-z) 


exp  i(  ?  )  -  exp  i(  -  ^  ) 


(20) 


(?F  ) 


1 


(^  ■  "zV  )  (^  ■1  TT  ) 


exp 


i^aT^  -  2TrfTjj 


1 


(f  -  ^) 


“T“ 

TT 


1 


(f  -  ^)  +  1 


TT 


+  0(-^) 

f 


TFT 


77 


4T 


exp 


[i( 


aT^  -  2TrfT]] 

~1 - 


4T 


+  0(4)  (21) 

f 
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7.  Sine-Squared  Pulse 


h(t)  =  sin^  (^  )  =  -  ^  ( 


.  TTt  .  TTt 

^  "T 

e  -  e  ) 


(22) 


G(f,  a)  ^  (-  -^  ) 


1 


1 


_  -  _ 

4(f-^)  4(f+^) 


+  exp 


-1 


1 


1 


2(f-^)  4(f-^-l)  4(f-^+l) 


4itT 


8.  Sine-of- 0-Squar ed  Puls e 


1 


?T3 

T 


exp 

j  {aT^  -  2^fT)] 

(f-^) 

+  0(4) 

f 


(23) 


h(t)  =  sin  0^(t);  0^(t)  =  it 


(T-t) 


2  1 


(24) 


G(f,  a)  s  (^  ) 


1 


1 


(f  +  4)  (f-  1)_^ 


4-nr  T 


exp  [i(aT^  -  ^T7fTj[' 


(f.^ 


+  (-  -4) 

Sir 


a  + 


a  - 


T 

VTJ 


(f+4) 


(f  -  y)^  j 


+  0(^) 


1 


(f 


+  (■ 


4ir  T 


?~J 


— 2“?  ^ 
4tt  T 


exp  -  2TTfTi 


+  (-  -4-) 

4-nr  T 


+  (-4-) 


TT  T 


i?47^ 

T 


"  'W’ 


T 


(25) 
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The  sine-of- 0~ squared  pulse  is  a  rather  interesting  one,  since  it 
resembles  certain  physically  realizable  pulses  having  rapid  rise  times  and 
well-defined  tails.  Moreover,  the  leading  edge  (the  part  to  the  left  of  the 
peak)  or  the  trailing  edge,  or  both,  could  be  fitted  together  with  other 
functions  to  approximate  other  pulses  having  one  or  both  of  the  properties 
mentioned  above.  Because  of  the  practical  potentialities  of  the  sine-of-0- 
squared  pulse,  a  detailed  version  of  this  pulse  is  presented  in  Fig,  2. 
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3-41-7664  1 


RECTANGULAR  PULSE 
h(t) 


PARABOLINEAR  PULSE 


- - 1 

- - - ii- 


T/2 


h(t)  =  1 


TRIANGULAR  PULSE 


h,(t)=  Y  in  0<  K ; 
h2(0  =  (Y)  (-f+T)  in|<KT. 


PARABOLIC  PULSE 
2+-^ - A 


h(t) 
1| — 


0  t:  T/2  \  T 


2 


Lh(t) 


T/2 


1  /  V  .  2  TTt 

h(0  -  sin  -y-. 


SINE-OF-e-SQUARED  PULSE 


hj(t)  =  in  0  <  K  =  --y--  ; 

,  h(t) 

-  1  in  <  K  ^  /* 

^  in  <  t  <  T 

1 

1 

0 

T/2  T 
h(t)  =  sin  0^(t)  ; 

02(f) 


Fig.  1 
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1,(1 

1  -  f 

■'•t 

t 
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tt/2 
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1 . 0000 

2tt/3 
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2 
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-f)^}  ; 

'■•T 

Fig.  2 


10 


C.  Description  of  the  Methods  Used  in  the  Evaluation  of  G(f,  of) 


The  evaluation  of  the  integral  for  G(f,  Qf)  in  Eq.  (3)  depends  primarily 


on 


related  integral 

t 

n 

^  exp  [ 

i  -  2TTft|J 

(26) 


s 

n 

subject  to  the  inequalities  stated  in  Eqs.  (6)  and  (7).  For  convenience, 
will  first  be  rewritten  in  another  form  by  completing  the  square  in  the  ex¬ 
ponent  and  then  introducing  the  substitutions 


^  =  A  -  ;  A  =  2fC  =  )f  ;  C  =  ;  a  =  ir/ZC^  ; 

s  Of  s 

V  =  A-2B  ;2B  -  ^;v  =  (2C)(f - 2-  )  (27) 

n  n  nC  n  tt" 

t  at 

u  =  A  -  2A  ;  2A  =  u  =  (2C)(f - ^  )  . 

n  n  nC  n'^'  tt  ' 


After  completing  the  procedure  just  described,  one  obtains,  in  place  of 
D^^(f,  0^),  the  following  integral  expression; 

i  TT  A^ 

^  ®  J^(A,  C)  ;  (28) 

where 

V  ^ 

n  .  TT  ^  ^ 

(  "  7  ^ 

J^(A.C)  =  J  e  de  . 

u 

n 

I  (and  hence  D  )  can  now  be  evaluated  provided  the  Fresnel-type  integral, 
on  on  ^  /  X  o 

J^,  in  the  above  equation  can  be  evaluated  subject  to  the  inequalities  stated 
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in  Eqs.  (6)  and  (7).  During  the  course  of  the  work  being  reported  here,  it 
was  found  that  one  could  evaluate  J^(A,  C)  by  employing  successive  integrations 
by  parts,  while  utilizing  certain  properties  of  the  curve  known  as  the  Cornu 
Spiral*  For  example,  the  first  two  integrations  by  parts  are  carried  out  in 
the  following  manner: 


V  V  ^ 

n  TT  n  ^  TT 

J^(A,  C)  =  ^  (iiT^e  ^  d^)  5  ^  d(e  ^  )  ; 


u 


n 


u 


n 


.  TT  ^2  .  TT 

-1  -1-2  ^7  ^  ^ 

w  =  (i'H'^)  ;  dw  =  “(i'H')  ^  d^  ;  dv  =  d(e  )  ;  v  =  e 


.  TT  Z  .  TT  2  n  1-^  5 

'  ITTV  ITTU  J  t 

n  n  (iTr)  § 

u 

n 


n 


(29) 


.  TT  ^2  .71^2 

-2  -3  -2  -4  ^7  ^  ^ 

w  =  (iTr)  ^  ;  dw  -  -(iTr)  (3^  )d^;  dv  =  d(e  )  ;  v  =  e 


J  (A,  C)  = 

n'  /  ITTV 


1 


TTTZ^ 

n  (ITT)  V 

n 


^  .  TT  2 

V 

2  n 

*  O  » 

^  KZ 

iTTu  /•  \2  3 

n  (iTT)  u 

/ 

2  n 


n 


y  - 


u 


(iir)  ^ 


•  ^  t2 
^2  ^ 

Continuing  in  this  manner,  each  time  setting  dv  equal  to  d(e  )  and  w 

s  t 

equal  to  the  rest  of  the  integrand,  one  obtains,  after  the  (M+1)  integration 
by  parts,  the  following  expansions  for  the  functions  and 


*See  Section  E. 
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where 


M 

J^(A,  C)  =  2 

m  =  0 


n 


mn  n 


.  .  TT  2  , 
l(  -X  V 

'2  n 


M 

I  (A,  C)  =  2 

m  =  0 


mn 


n 


(■ 


+  C  exp  ( -  i  ^  A^  )  R 


Mn 


(30) 


The  functions  and  in  Eqs.  (28)  and  (30)  can  now  be  used  in 

the  evaluation  of  certain  important  related  integrals,  namely: 


t 

n 

cS  =  1  a  exp  fifcrt^  -  27rft]1dt  =  a  D  (31) 

on  j  on  ^  L  l  JJ  on  on 

s 

n 

*The  convergence  properties  of  this  asymptotic  expansion  are  discussed  in 
the  appendix. 
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(31) 


a 

In 


2n 


n 

^  aint  exp  [i  -  2TTft}]dt  = 


dD 


on 


“ST 


n 


n 


.1' 


t  exp 


i[at^  -  2Trft]]dt  =  a2^(^) 


.  a^D 

^  on 


af 


where  a  ,  a,  ,  and  a^  are  constants.  One  or  more  of  these  integrals 
occurs  in  the  evaluation  of  the  spectrum  of  each  of  the  eight  pulses  previously 
listed.  For  example,  all  three  are  applicable  when  the  pulse  shape  is  either 
parabolinear  or  parabolic,  while  only  the  first  two  apply  in  the  case  of  the 
triangular  or  the  trapezoidal  pulse.  The  above  three  integrals  can  be  ex¬ 
pressed  in  terms  of  I  and  the  first  two  derivatives  of  I  with  respect  to 
A  with  the  aid  of  the  following  relationships: 


A  =  2£C  ;  D  (£,  O')  =  D  — I_)=I  (A,  C)  ; 

on'  '  on'ZC 

ao  91  Q  A  91  a^D  .  a^i 

on  _  on  oA  _  on  on  _  on 

-TJ—  =  “3X"  ^  ^ 


(32) 


Equation  (31)  can  now  be  re-written  as  follows; 


ai 


a  =a  I  ;  a,  =a.  ;  a  (i£) 

on  on  on  In  in  tt  oA  Zn  Zn  tt 


P  a^i 
^  on 


aA 


(33) 


When  the  differentiations  indicated  in  Eq.  (33)  have  been  performed,  using 
the  expression  for  given  in  Eq.  (28),  one  obtains 
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-i-j-  A 


■^1  ^ 


^  =  (a  C)  fj  1 

on  '  on  '  L  nj 

=  (a,  C^)  Taj  +(-J-)J'1e  "  (34) 

In  '  in  L  n  '  ITT  '  nj  '  ' 

=  (a.  C^)  [(-  JL)J  +A^J  +  (-  ^)  A  J'  +  j"] 

2n  2n  L  i'”’  n  n  iir '  n  iir  nj 


-ij  A 


where 


aj 


n 


a^j 


Tl  -  _ tL  •  T  "  - 


a  A' 


and 


-iZa^  aj  -iZa^  iT] 

J*  e  =  --r;  e 


^  TT 


n  n 

=  e  -  e 


J"  e 
n 


-i-^-A 


r,  2  ^  .  TT  .2  .  . 

a  j  - 1  ^  A  1/^  IT] 

n  2  ...  n  ...  n 

- y-  e  =  (iTTv  )  e  -  (iTTu  )  e 

aA^  ^  ^ 


(35) 


A  ,  TT  .  2  2  ■’’■  /  2  2. 

A  =  V  +  -7^  £u  +-7>^;t(v  -  A)e.^;^(u  -  A  )  -  rj 

n  ^  n 


n  ^  -C  =  %  ^ -C-  ’  I  '"n 


n 


The  next  step  is  to  replace  in  Eq.  (34)  by  the  expansion  given  in  Eq.  (30) 

and  J‘  and  J’*  by  the  expressions  given  in  Eq.  (35).  The  result  of  making 

s 

these  substitutions  and  combining  terms,  while  employing  A  =  v^  4  in 

conjunction  with  terms  involving  v^  and  and  A  =  u^  4  in  conjunction 
with  terms  involving  u^  and  ,  is  shown  in  the  following  equations: 


5. 


on 


iT? 


(a  C)  V  e 
on  mn 


-  (a  C)  U  e 
on  mn 


'’]}  "  f’^on  } 


(36) 
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a 


In 


Mr  i7] 

=  <'  2  (a,  s  C)V  e  ^  -  (a,  t  C)  U  e  ' 

In  n  mn  In  n  mn 


(m  =  0  . 


M 
2 

m  =  l  L 


P  2 

(a,  C  )v  V  e  -  (a,  C  )u  U  e 
In  n  mn  In  n  mn 


WrinU37) 


a 


2n 


M  r  2 

2  (a,  s  C)V  e  "  -  (a,  t  C)U  e  " 

„  Wn  n  mn  tin  n  mn 

m  =  U  U 


+  J  2 

m  =  l 


Mr  2  2 

{2a,  s  C  )v  V  e  "  -  (2a  ,  t  C  )u  U  e 
Zn  n  n  mn  Zn  n  n  mn 


+ 


^  M  2  7 

V  /  3m-2  ^ 
m  =  Z 


'2  9  '2  9 

.  3.  Z  n  .  3v  Z  n 

(ao  C  )v  V  e  -  (a-5  C  )u  U  e 
Zn  n  mn  Zn  n  mn 


a^  C  irj 

+  {  ( - - )(V  e  -  u  e  " 

'  '  ITT  Mn  Mn 


('a„] 


(38) 


The  remainder  terms  r  ,  r,  ,  and  r^  in  the  above  equations  are  defined 

on  In  Zn 

as  follows: 


r  =  (a  C  e 
on  on 


■’■z  ^ 


i 

’■in  =  ("!„=" 


3  1  2 

r,  =  (a,  C  e  ){-^  +  A  ){R,  ,  )  ;  R,  , 

2n  2n  itr  Mn'  Mn 


V 


^  b 


•  ^ 

M+1  ,, 

2m+3  ^  ’ 


U 


n 


’^M+l 


{-1){1){3) 


{iT7) 


-  {2M-1){2M+1) 
M+2 


(39) 


V 

n 


as  t  n't 

{2C){£  -  -^);  =  A  -  -^  =  {2C){£  -  -^  ) 
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The  functions  o/  ,  ,  ,  and  in  Eqs.  (36),  (37),  and  (38)  can  now 

on  In  cn  \  \  / 

be  combined  in  order  to  permit  the  computation  of  the  FM  spectra  of  three 

basic  waveforms,  each  one  defined  in  the  interval  s  <  t  <  t  .  The  three 

n  n 

waveforms  referred  to  above  correspond  to  the  following  definitions  of  the 
function  g(t)  appearing  in  Eq.  (1): 


g(t)  -  h^^(t) 


-  a 

on 


constant  in  s  <  t  <  t  (Rectangular  Pulse 
n  n  '  ^  .. 

oegment) 


g(t)  h^(t) 


a 

on 


in  s  <  t  <  t  (Linear  Pulse 
n  n  c 

b  egment) 


(40) 


g(t)  =  h  (t)  =  a  +  a,  t  +  t  in  s  <  t  <  t  (Parabolic 
^  n  on  In  Zn  n  ^  ^  „  . 

Pulse  oegment) 


3-41-7666 


h  (t)  =  a 

n  On 


RECTANGULAR 

SEGMENT 


h  (t)  =  a  4-  a  t 

n  On  1n 


LINEAR 

SEGMENT 


h(t)=a  +at-i-a  t 

n  On  In  2n 


CASE  HI 

PARABOLIC 

SEGMENT 


Fig.  3 


*See  Fig.  3  and  Eqs.  (40)  through  (43). 
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The  results  can  be  summarized  as  follows: 

Case  I  Rectangular  Pulse  Segment  h^(t)  =  a^^  in  ^  <  t^ 


a 

on 


M 

E 

m=0 


fh  (s  ) 

I  n  n 


CV 


mn 


-  h  (t  ) 
n  n 


CU 

mn 


(41) 


wher  e 


h  (s  )=h  (t)=a 
n  n  n  n  on 


Case  II  Linear  Pulse  Segment 


a  +aitins  <t<t 
on  In  n  n 


(41a) 


M  , 

^  -  E  (s  ) 

on  In  ^ 1  n  n 

m=0 


CV  e 
mn 


M  r 
+  E 
m  =  l 


h*  (s  ) 
n  n 


C  V  V  e 
n  mn 


n 

i- 


h  (t  ) 
n  n 


CU  e 
mn 


h*  (t  ) 
n  n 


■7  IT]  ■ 

TT  ^ 

C  u  U  e 
n  mn 


+  +  r . 

on  1  n 


(42) 


where 


h  (s  )  =  a  +  a,  s  ;  h  (t  )  -  a  +  “  ^In 

n  n  on  Inn  n  n  on  Inn  s  n  in 


Case  III  Parabolic  Pulse  Segment  h^(t)  =  a^^  +  a^^t  +  a^^t  in  <  t  <  t^ 


^  +  ^1  +  ^  ^  [h  (s  ) 

on  In  2n  n  ^  ^ 

m  =  U  L 


CV  e 
mn 


h  (t  ) 
n  n 


117 


CU  e 
mn 


n 


M  f 
+  2 
m  =  1 


h'  (s  ) 
n  n 


■7  1 

^  2  ,  ^  n 

C  V  V  e 
n  mn 


-  h*  (t  ) 
n  n 


^  2  n 

C  u  U  e 
n  mn 


*07 

3  2  n 

C  V  V  e 
n  mn 


-  h**(t  ) 
n  n 


C^u^U  e 

n  mn  . 


h"(s  ) 
n  n 

~ZH“ 


i^  1 

n 


h"(t  ) 

n  n 

ZttI 


Q  ‘ 

/-^TT  “ 

c  u,  ,  e 

Mn 


'+  ''on  +  'in  +  '2„ 


] 


(43) 
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where 


2  2 
h  (s  )  =  a  +  a,  s  +  a^  s  ;  h  (t  )  =  a  +  a,  t  +  a^  t 
n  n  on  Inn  ^n  n  n  n  on  In  n  2n 


(43a) 


h*  (s  )  =  a,  +  2a^  s  ;  h'  (t  )  =  a,  +  t  ;  h"(a  )  =  h"(t  )  =  2a^ 
n  n  In  2n  n  n  n  In  2n  n  n  n  n  n  2n 
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D.  Outline  of  the  Problem  of  Evaluating  the  FM  Spectrum,  G(f,  Q^),  of  a 
Pulse  Describable  as  a  Piecewise  Linear  Function  of  Time  in  the 
Interval  0  ^  t  ^  T 

To  illustrate  the  use  of  the  formulas  in  Eqs.  (41)  through  (43a)  in 
calculating  the  FM  spectrum  of  a  useful  function,  the  problem  of  calculating 
the  spectrum  of  a  pulse  which  can  be  approximated  by  a  piecewise  linear 
function  will  now  be  considered  briefly.  In  this  problem,  using  the  definition 
of  h(t)  for  the  piecewise  linear  case,  as  given  in  Eqs.  (8)  and  (9),  one  is 
faced  with  the  evaluation  of  the  following  integral: 

Git.c)  .  yh(t)  ; 

o 

h(t)  h  (t)  in  (n-l)(^)  ^  t  ;  h(t)  =  a  +  a,  t 

'  n  N  In  ’  '  '  on  In 


which  is  a  special  case  of  the  general  integral  defined  in  Eqs.  (1),  (2),  and 
(3)  of  this  report.  In  performing  the  required  evaluation  of  the  integral  in 
Eq.  (44),  Case  11,  as  described  in  Eqs.  (42)  and  (42a),  is  the  case  which 
applies.  Equation  (44)  can  thus  be  rewritten  in  the  following  manner; 


G(f.c^) 


N 

Z  - 
n=l 


dt  + 


S 

n 


,  2  {2  . 

n  =  l  ^ 


(45) 


wher  e 


s 

n 


(n-l)(^)  = 


nT 

N 


n( 


nT 


A  formal  expression  for  the  spectrum  described  in  Eq.  (44)  can  now  be 
obtained  by  substituting  in  Eq.  (45)  the  asymptotic  expansion  of  the  function 
‘^on  ^n  given  in  Eqs.  (42)  and  (42a),  Once  again,  subject  to  the  assumption 
stated  in  Eq.  (7),  a  suitable  approximation  for  G(f,  Q^)  can  be  obtained  by 
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retaining  only  the  first  non- vanishing  term  of  each  of  the  expansions  given 
in  Eq.  (42).  The  addition  of  terms  such  as  these  for  all  of  the  N  intervals 
into  which  the  interval  0  ^  t  ^  T  has  been  divided  of  course  entails  additional 
mathematical  complications.  However,  these  are  not  too  severe,  especially, 
as  is  often  the  case,  when  the  pulse  shape  h(t)  can  be  adequately  approximated 
by  a  piecewise  linear  function  for  which  the  number,  N,  of  sub-intervals  is 
between  5  and  10.  For  values  of  N  smaller  than  about  5,  the  computations 
required  are  relatively  simple.  In  particular,  Case  II  was  employed  with 
a  minimum  of  difficulty  in  the  calculation  of  the  spectrum  of  the  triangular 
pulse,  the  results  of  the  calculation  having  been  presented  earlier  in 
Eq.  (13).  Inthiscase  only  two  sub-intervals  were  used  (i.  e.  ,  N  =  2). 

Similarly,  the  spectrum  of  the  rectangular  pulse,  as  given  in  Eq.  (11),  was 
calculated  with  the  aid  of  the  result  given  for  Case  I,  with  N  =  1  ,  while 
the  spectra  of  the  trapezoidal,  parabolic,  and  parabolinear  pulses,  as 
given  in  Eqs.  (17),  (15),  and  (18),  respectively,  were  calculated  using  the 
appropriate  combinations  of  Cases  I,  II,  and  III  and  at  most  a  value  of  N  -  3. 
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Geometrical  Interpretation  of  the  First  Term  of  the  Expansion  of  J  (A,  C) 


n 


It  was  stated  earlier  that  certain  properties  of  the  curve  known  as 
the  Cornu  Spiral  could  be  utilized  in  conjunction  with  the  integration-by-parts 
procedure.  While  it  is  not  at  all  necessary  to  invoke  the  geometrical 
characteristics  of  the  Cornu  Spiral  in  this  connection,  nevertheless  by  so 
doing  one  can  provide  convenient  geometrical  interpretations  of  the  various 
terms  in  the  expansion  of  J^(A,  C)  given  in  Eq.  (30),  and  thus  provide  a 
much  more  satisfying  analysis  than  one  based  on  the  purely  formal  integration- 
by-parts  procedure  outlined  above.  In  the  following,  only  the  first  term  in 
the  expansion  and  its  geometrical  interpretation  will  be  discussed. 

The  coordinates,  x  and  y,  the  slope,  ,  the  second  derivative,  — -y  > 
and  the  radius  of  curvature,  p,  of  the  Cornu  Spiral  at  a  point  P(x,  y)  are 
given  in  the  following  equation. 

u  u 

x(u)  =  ^  cos  ~  d^  ;  y(u)  =  ^  sin  ^  ^  d^  ;  u  =  OP 


J(u)  =  X  +  iy 


.  7T  ^  Z 

e  d^  =  OP  (see  Fig.  4) 


dx  =  cos  ^  u  du  ;  dy  =  sin  ^  u^  du 


(46) 


dy 


=  y*  =  tan 


77 


du 


sec  u 


d^y  „  _  ,  dy' 

- T  =  V  =  ^ 

1  ^  du 

dx 


77  u  sec 


3  77  Z 


Slope  at  P(x,  y)  ; 


p  =  Radius  of  Curvature  at  P(x,  y) ;  p  =  1  +  y* 


3/2  ,,  ,  _1 

^  TTU 


The  function  J,  as  shown  in  Fig.  4,  is  the  vector  OP  from  the  origin  0  to 
the  arbitrary  point  P  on  the  spiral.  Similarly,  the  function  in  Eq.  (Z8) 
is  the  vector  shown  in  Fig.  4,  or  equivalently,  the  vector  difference 


ZZ 


X  =  x(u)  =  ^  cos  ^  ^  ;  y  =  y(u)  =  ^  sin  ^  ; 

—  °  d  nu^ 

u  =  OP  =  arc  length  at  P  ;  ^  =  tan  =  slope  at  P  ; 

J(u)  =  X  +  iy  =  C  exp  {i  ^  d^  =  OP  ■ 


C  =  center  of  curvature  of  spiral  at  P 


p  =  CP  =  —  =  radius  of  curvature  at  P 

TTU 


Fig.  4 
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3-41-7668 


.  i  tt/2  a f. 

J  (A,C)  =  \  e  ' 
n  Ju^ 

^  r"e'"/2  5 


=  OQ  -OP  H 

n  n 


P  Q 

n  n 


Fig.  5 
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Fig.  6 
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PP^~PC  ;  QP^~QD  ; 
nn  n°o  nn 


J  (A,C)  E  P  P^  -  Q  P^ 

n'  n  n  00 


P  C  -  Q  D 
n  n  n  n 


Fig.  7 
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- ^  ^ 

-  ^n^oo  where  P^=  ^oo^’  limiting  point, 

the  point  at  which  the  arc  length,  measured  from  0,  is  infinite.  It  is  clear, 
after  knowing  the  meaning  of  C)  and  noting  the  definitions  of  the  upper 

and  lower  limits,  v  and  u  ,  respectively,  given  in  Eq.  (27),  that  when  the 
inequality  in  Eq.  (7)  is  satisfied,  both  ends  of  the  vector  (i*  >  the 

function  J^(A,  C))  lie  far  within  the  wound-up  portion  of  the  spiral.  In  this 
portion  of  the  spiral,  the  center  of  curvature,  C,  of  the  spiral  and  the 
limiting  point,  sire  very  nearly  coincident,  so  nearly  so  in  fact,  that  in 

a  great  many  applications  (but  not  in  the  present  one)  the  following  is  an 
excellent  approximation  for  the  integral  J^(A,C).  (See  also  Fig.  7.) 

n  .  IT 

(■  ^7  ^  ^ 

J  (A,  C)  =  \  e  d|  =  P  P  -  Q  P 

n'  '  J  ^  n  oo  n  oo 

u 

n 


PC  -  Q  D  = 
n  n  n  n 


.(  TT 

T7  2) 

exp 

7''niJ 

_ 

TTU 


TTV 


n 


(47) 


where 


TTUn 


and 


-  are  the  radii  of  curvature  at  P  and  Q  respectively. 

TTV  n  n  ^  ' 


Equation  (47)  was  obtained  by  first  noting  that  the  lengths  of  the  vectors 


P  C  and  Q  D  are  simply  the  radii  of  curvature  — ^ 
n  n  n  n  ^  '  ttu 


and 


1 


TTV^ 


,  respectively. 


n  n 

at  the  two  points  P  and  Q  .  Next,  it  is  noted  that  the  phase  angles  of  the 

n  -J  n 

IT  TT  d  .  TT  ,  TT  ^ 

two  vectors  are  '2  “Z and  2  "Z ’  respectively.  For  example,  the 
IT  Z  ^ 


term  -j- u  is  simply  the  angle  of  inclination  of  the  spiral  at  P  ,  as  given  in 
general  in  Eq.  (46).  The  remaining  angle,  >  appears  because  the  spiral  is 
being  approximated  at  P^  by  a  circle  centered  at  C^,  and  thus  P^C^  is 
perpendicular  to  the  spiral.  The  above  ideas  are  illustrated  schematically 
in  Fig.  7  . 

Upon  comparing  Eq.  (47)  with  the  third  line  of  Eq.  (29),  one  sees  that 
the  approximation  given  in  Eq.  (47)  is  none  other  than  the  FIRST  TERM  in 
the  expansion  J^(A,C).  Now,  the  locus  of  the  center  of  curvature  of  the 
spiral  is  called  the  EVOLUTE  OF  THE  SPIRAL.  One  can  easily  show  that 
the  second  term  in  the  expansion  (terms  in  and  — j)  is  related  to  a 


u 


n 


n 
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pair  of  vectors  pointing  from  to  the  center  of  curvature  of  the  EVOLUTE 

at  C  and  from  D  to  the  center  of  curvature  of  the  EVOLUTE  at  D  . 
n  n  n 

Continuing  in  this  manner,  one  can  give  geometrical  interpretations  to  all 
of  the  remaining  terms,  although  the  geometry  involved  becomes  increasingly 
complicated  with  each  succeeding  term. 
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APPENDIX  A 


Some  Aspects  of  the  Convergence  Properties  of  the  Asymptotic  Expansions 
of  the  Fresnel  Integrals  J^(A,C)  and 

In  Section  C  an  integration-by-parts  procedure  was  used  to  obtain  a 
purely  formal  asymptotic  expansion  of  the  integrals  J^(A,C)  and  I^^(A,C) 
defined  in  Eq.  (28).  The  results  obtained  after  one  and  two  integrations 
were  given  in  Eq.  (29).  The  general  result  of  employing  M+1  (recall  the 
index  m  runs  from  0  to  M)  integrations  was  shown  in  Eqs.  (30)  and  (30a). 
However,  nothing  was  said  at  the  time  concerning  the  convergence  properties 
of  the  expansions  obtained  in  this  purely  formal  manner.  In  particular, 
nothing  was  said  concerning  the  magnitude  of  the  error  incurred  by  neglecting 
the  remainder  term  Eq.  (30)  and  using  the  M+1  terms  in  the  ex¬ 

pansion  to  approximate  the  integral. 

The  purpose  of  this  appendix  is  to  obtain  an  upper  bound  on  the  magni¬ 
tude  of  the  error  introduced  by  discarding  all  but  a  finite  number  of  terms 
in  the  asymptotic  expansion  of  the  function  (and  hence  given  in 

Eq.  (30).  The  decidedly  more  difficult  problem  of  estimating  a  lower  bound 
on  this  error  will  not  be  treated  here,  but  it  certainly  merits  attention  at 
some  later  date.  To  aid  in  the  estimation  of  the  upper  bound,  the  expression 
for  J^(A,C)  given  in  Eq.  (30)  will 

V  ^ 

n  .  TT 

C  ^7  ^ 

J^(A,C)  =  ^  e  de 

u 

n 


now  be  rewritten  in  the  following  form: 
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.  7T  Z  .  7T  2 
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(2M+l)b 
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ITT 
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=  (2C)(f  -  ^);  =  (2C)(f 


(A -2) 


In  what  follows  use  will  also  be  made  of  the  following  important  inequalities: 

1 

2M  +  1 


u 


n 


V 


which  follow  immediately  from  the  definitions  of  u  and  v  given  above  and 


n 


in  Eq.  (27)  and  from  the  way  in  which  the  limits  of  integration  s^  and  t^ 
appear  in  the  definition  of  the  basic  integral  given  earlier  in 

Eq.  (26): 


n 

D^^(f ,  O')  =  ^  exp  -  27Tft|j  dt 


(26) 


.  St 


A  comparison  will  now  be  made  of  the  magnitude  of  the  (M  +  1)  term 

in  the  asymptotic  expansion  of  J^,  and  the  magnitude  of  the  remainder 

integral  one  obtains  after  performing  M+1  integrations  by  parts.  The  term 


can  be  rewritten  as  follows: 
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The  magnitude,  that  is,  the  absolute  value,  of  this  complex  vector  function 
IS  seen  immediately  to  be: 
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where  use  was  made  of  the  inequalities  stated  in  Eq.  (A-3).  Required  also 
are  the  maximum  and  minimum  values  of  the  function  |  .  They  are: 
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Mn‘  MAX 


TMTl 
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“ZMTT  ' 

V 
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(A- 6) 


obtained  by  setting  i  (v  -  u  )  in  Eq.  (A-5)  equal  to  (2j  +  l)7r  (j  =  0,  ±1  ,  ±2,  .  .  ) 

^  n  n 

and 
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Mn '  MIN 
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obtained  by  setting  ^  (v  -  u  )  -  Zj-n  (j  -  0,  il,  ±2,  .  .  .  ). 

^  n  n 

The  magnitude  of  the  remainder  integral  can  best  be  examined 

by  rewriting  this  integral,  as  defined  in  Eqs.  (A-2)  and  (30a),  in  the  following 
form: 


^Mn 
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Because  of  the  difficulties  involved  in  integrating  the  above  expression,  a 

closed-form  expression  for  the  magnitude  of  unattainable,  although 

approximate  integration  techniques  can  be  invoked  if  necessary  to  give  a 

result  close  to  the  true  value.  No  attempt  in  this  direction  will  be  made 

here,  however.  Instead,  an  estimate  of  an  upper  bound  (that  is,  a  maximum 

value)  on  the  magnitude  of  will  be  deemed  sufficient  to  meet  the  needs 

of  the  present  analysis.  This  upper  bound  can  be  easily  obtained  if  it  is 

observed  that  for  the  range  of  values  of  u  and  v  under  consideration  (as 

determined  bv  the  definitions  of  u  and  v  in  Eq.  (27)  and  by  the  fundamental 
^  n  n  ^  ^ 

inequalities  stated  in  Eqs.  (6)  and  (7))  the  phase  of  the  complex  exponential 

.  TT  uZ 

term  e  in  Eq.  (A-8)  is  an  extremely  rapidly  oscillating  function  of  ^ 

and  thus  in  the  region  of  integration,  u^  <  ^  <  v^,  there  occurs  a  considerable 


n 


amount  of  destructive  interference  between  the  incremental  complex  vectors 
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^2M+2  ® 


di 


involved  in  the  integration  process. 


On  the 


basis  of  the  above  ideas,  one  can  immediately  deduce  the  following  very 
strong  inequality: 
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Combining  the  results  presented  in  Eqs.  (A-6),  (A-7),  and  (A-9),  one 
now  obtains  the  following  relationship  between  the  various  magnitudes: 


«  ( 


~ZM+1 

u 

n 


V 


M' 

■zM+r 


)  = 


T  < 

Mn '  MIN 


Mn'MAX 


(A-10) 


From  this  equation  one  obtains  immediately  the  required  relationship 

defining  an  upper  bound  on  the  magnitude  of  the  error  introduced  by  neglecting 

the  remainder  term  appearing  in  the  asymptotic  expansion  of  the 

function  J  ,  namely: 

n  ^ 
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Mn  • 
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Mn ' 


(A-ll) 


Equation  (A-ll)  shows  that  the  magnitude  of  the  remainder  term  is 

s  t 

always  much  less  than  the  magnitude  of  the  (M+1)  term  in  the  series,  and 
a  closer  examination  reveals  that,  for  a  fixed  value  of  M,  the  magnitude  of 
the  remainder  (and  hence  of  the  error)  grows  continually  smaller  with 
increasing  u^  and  v^.  It  is  evident,  then,  that,  by  terminating  the  series 
at  a  term  whose  absolute  value  corresponds  nearly  to  the  maximum  error 
tolerable  in  the  calculation  of  the  function  J^,  one  can  establish  an  upper 
bound  on  the  error  incurred  by  neglecting  the  rest  of  the  terms.  The  same 
considerations  of  course  also  apply  to  the  calculation  of  the  related  function 
I^^  in  Eq.  (30).  Finally,  these  ideas  can  be  extended  relatively  easily  to 
the  discussion  of  the  convergence  properties  of  more  general  asymptotic 


32 


expansions  such  as  those  presented  in  Eqs.  (4Z)  and  (43),  and,  in 
expansions  of  integrals  of  the  following  general  type: 

G  (f,  o^)  =  r  (a  +  a,  t  +  a^  t^ )  e^^  ^TrftJ  ^ 
n  J  '  on  In  Zn 

s 
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